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ABSTRACT 

We study the invariant properties of the convection equation du/dt + a du/dx = 0 with respect to 
spatial reflection, time reversal, and space-time inversion. Generally, a finite-difference analogue 
of this equation may possess some or none of these properties. It is shown that, under certain 
conditions, the von Neumann amplification factor of an analogue satisfies a special relation for 
each invariant property this analogue possesses. Particularly, an analogue is neutrally stable and 
thus free of numerical diffusion if it possesses the invariant property related to space-time 
inversion. It is also explained why generally (i) an upwind scheme possesses neither the invariant 
property related to spatial reflection nor that related to space-time inversion, and ( ii ) an explicit 
scheme possesses neither the invariant property related to time reversal nor that related to space- 
time inversion. Extension to the viscous case and a remarkable connection between the current 
work and a new numerical framework for solving conservation laws are also discussed. 


1. INTRODUCTION 


Many physical equations are invariant (i.e., they do not change their forms) under certain 
transformations. As a simple example, consider the convection equation 


du du n 

dt ox 

where a is an arbitrary real constant. Eq. (1 . 1) is equivalent to 

du , du 

a7 + 3 a7 = 0 


if 


or 


(0 

(«) 


x' = -x , t' - t , and a' = -a 
x' = x , t' = -t , and a' = -a 


(U) 

( 1 . 2 ) 

(1.3) 

0.4) 


In other words, Eq. (1.1) is mapped into the same equation under the mapping 

x — > x' , t — » t' , and a — » a (1.5) 

if x', t ' , and a' are defined either by Eq. (1.3) or Eq. (1.4). Since the mapping defined by Eq. 
(1.5) with 

(Hi) x' = -x , t' = -t , and a’ = a (1.6) 

can be considered as the product of the first two mappings (in either order), Eq. (1.1) is also 
mapped into itself under the third mapping. 
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Since the mapping (x,t) -> represents (i) spatial reflection with respect to the reference 

plane x = 0 if x = -x and t' = t, (ii) time reversal with respect to the reference plane t = 0 if x' = 
x and t' = -t, and (iii) space-time inversion with respect to the origin if x = —x and t' = -t, Eq. 

(1.1) is said to possess invariant properties with respect to spatial reflection, time reversal, and 
space-time inversion respectively. 

As a result of the above properties, solutions to Eq. (1.1) also possess similar invariant properties. 
Let F be a function of a real variable. Then, for any a, 

u = u Q (x,r,a ) d - F(x-at) (1.7) 

is a solution to Eq. (1.1) if the derivative of F exists. Since (i) u = u 0 (x,f,a ) is a solution to Eq. 

(1.1) if and only if u = u 0 (x' is a solution to Eq. (1.2), and (ii) Eq. (1.1) is equivalent to Eq. 

(1.2) if x', t', and a' are related to x, t and a by either Eq. (1.3) or Eq. (1.4) or Eq. (1.6), one 
concludes that 

u = u 0 (-x,t;-a) , (1.8) 

u = u 0 (x,-t;-a ) , (1.9) 

and 

u = u 0 (-x,-t;a) (1.10) 

are also solutions to Eq. (1.1) (Note: do not alter the sign of a in Eq. (1.1)). The above 
conclusion can also be verified directly by using the facts that u 0 (-x,t\-a) = F(-(x-at)), 
u 0 (x, -t;-a) = F(x—at), and u 0 (-x, - t;a ) = F(-(x-at)). 

Numerical schemes generally are constructed without considering the invariant properties of the 
physical equations to be solved [1]. As a result, numerical solutions generally do not share the 
same invariant properties of physical solutions. Using numerical analogues of Eq. (1.1) as 
examples, it will be explained in this paper how violation of the invariant property with respect to 
space-time inversion is related to the numerical diffusion of these analogues. Particularly, it will 
be shown that a two-level constant-coefficient difference analogue of Eq. (1.1) is neutrally stable 
and thus free of numerical diffusion if it preserves the invariant property related to space-time 
inversion. 

The remainder of this paper is briefly described as follows: In Section 2, we investigate the 
invariant properties of several well known numerical analogues of Eq. (1.1). It is shown that the 
Wendroff scheme Eq. (2.1) and the Crank-Nicolson scheme Eq. (2.11) are invariant under (i) 
spatial reflection, (ii) time reversal, and (iii) space-time inversion, in a sense to be defined. Note 
that invariance under any two of (i) — (iii) implies invariance under the third. We also show that 
the invariance of a scheme under any one of (i) - ( iii) cannot occur unless the configuration of its 
stencil satisfies a certain necessary condition. This explains why generally an upwind scheme 
cannot be invariant under either spatial reflection or space-time inversion while a two-level 
explicit scheme cannot be invariant under either time reversal or space-time inversion. 

In Section 3, we consider an arbitrary two-level constant-coefficient finite-difference analogue of 
Eq. (1.1). It is shown that the von Neumann amplification factor of this analogue must satisfy 
certain relation if it is invariant under any one of the transformations referred to above. As a 
result, one can show that a scheme must be neutrally stable if it is invariant under space-time 
inversion. 

In Section 4, the discussion is extended to the convection-diffusion equation Eq. (4.1). We also 
consider the Leapfrog/DuFort-Frankel scheme Eq. (4.22), which involves three time levels and 
also is invariant under all three transformations referred to above. 
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In Section 5, we summarize the key results obtained in the current study and explain how they 
can be used in a more complicated situation. We also show that, contrary to traditional two-level 
explicit scheme, a new two-level explicit analogue of Eq. (4.1), which is constructed by using a 
numerical frameworic currendy being developed by Chang and To [2], is invariant under spatial 
reflection, time reversal and space-time inversion. It is also shown that there is a remarkable 
similarity between the forms of the amplification factors of the new scheme and the 
Leapfrog/DuFort-Frankel scheme. The implication of this similarity is also discussed. 


2. INVARIANT PROPERTIES OF NUMERICAL SCHEMES 


A two-level difference scheme which can be used to solve Eq. (1.1) is the Wendroff scheme 
[p.503, 3], i.e., 

(l+v)(u$ -up + ( 1 - v)(«7- 1 -u] +l ) = 0 (2.1) 


where 


def aAt 

V 

Ax 

is the Courant number. We assume that Eq. (2.1) is valid for j,n = 0, ±1, ±2, 


( 2 . 2 ) 


Given any pair of integers j 0 and n 0 , there is one equation in Eq. (2.1) with j =j 0 and n =n Q . 
The image of this equation under the mapping SR defined by m " -» ulj (J, n =0, ±1, ±2, • • • ) 
and v -> -v is 

(l-v)( M ^ + 0 , + D-«-; 0 ) + (I+vXk^ 1 - u -Uo+l) ) = 0 (2-3) 

Eq. (2.3) can be rewritten as 

( 1 + v)(«!!^ 1 + i) + i -u^, 0+1 )) + (1 - v)(m"q^ 1 +1) — «"^ 0+1)+1 ) = 0 (2.4) 

A comparison between Eq. (2.4) and (2.1) reveals that the image of the original equation is also 
oneofEq. (2.1) with j =—(Jo+l), n=n q and the same Courant number v. 


Since j 0 and n 0 are arbitrary integers, one concludes that the system of equations represented by 
Eq. (2.1) is mapped into itself under the mapping SR. In other words, the system of equations 
represented by 

(1 -v)(M^Vi) -M-j) + (1 +v)(m _/ 1 -M-o+ 1 )) = 0 , j,n=0, ±1,±2, ••• (2.5) 

is identical to that represented by Eq. (2.1). In this paper, a scheme with this property is said to 
be invariant under spatial reflection. 


Let R\ be a set of real numbers such that the negative of an element is also an element. For any 
ve R j , let 

u] = u 0 (j,n;v ) , j,n = 0, ±1, ±2, ■ - • (2.6) 

be a solution to Eq. (2.1). Since Eqs. (2.1) and (2.5) are identical, Eq. (2.6) is also a solution to 

Eq. (2.5), i.e., 

( 1 - v) [ «o(“0+l).n+Lv) — uo(—j,n,v) ] 

+ ( 1 + V ) [ u 0 (-j,n+l ;v) - «o(-0' +l),rt ;v) ] = 0 , j,n =0, ±1, ±2, ••• (2.7) 

Since Eq. (2.7) is valid for any ve/? j , it is also valid if v is replaced by -v. A comparison 

between this new form and Eq. (2.1) reveals that, for any ve R \ , 

Uj = «o(-;.«;-v) , j,n =0, ±1, ±2, ••• (2.8) 
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is also a solution to Eq. (2.1). Obviously, this property is shared by any scheme which is 
invariant under spatial reflection. Since j is the numerical analogy of x and v = aAt/Ax, this 
property is similar to the property that Eq. (1.8) must satisfy Eq. (1.1) if Eq. (1.7) does. 

Similaiiy, it can be shown that the system of equations represented by Eq. (2.1) is mapped into 
itself under the mapping TR defined by u" — > uj n (j,rt =0, ±1, ±2, • • • ) and v — » -v. In this 
paper, a scheme with this property is said to be invariant under time reversal. Also it may be 
shown that, for any veR i , 

u J = uoU, -« ;-v) , j,n = 0, ±1, ±2, • • • (2.9) 

must satisfy a scheme which is invariant under time reversal if Eq. (2.6) satisfies this scheme for 
any ve/?!. 

Let the mapping STI be defined by u* — > uJJ (J,n = 0, ±1, ±2, • • • ) and v — > v. Since any one 
of the three mappings SR, TR, and STI is the product of the other two mappings (in either order), 
the system of equations represented by Eq. (2.1) is also mapped into itself under the mapping 
STI. In this paper, a scheme with this property is said to be invariant under space-time inversion. 
Obviously, for any ve R p, 

Uj = Mo (-/,-« ;v) , j,n= 0, ±1,±2, ••• (2.10) 

must satisfy a scheme which is invariant under space-time inversion if Eq. (2.6) satisfies this 
scheme for any ve R \ . 

To proceed further, we now establish a necessary condition for a scheme to be invariant under 
spatial reflection. Consider the equation in Eq. (2.1) with j=jo and n = n 0 . Its stencil S 0 is 
formed by the mesh points A = (Jo> n oX B = 0'o+l’ rt o)> C = 0‘o+l.no + l)* a °d D = (j 0 ,n 0 +\) (see 
Fig. 1). The image of the above equation under the mapping SR has a stencil So formed by the 
mesh points A' = (-j o,«o). B' - (-Oo+l). n o). C' — (-Oo + l).«o+l)> and£>' = (-Vo,«o+l)- 

Since the Wendroff scheme is invariant under spatial reflection, the image referred to above is 
also one of Eq. (2.1) (with j = -O’o+l) and n = n 0 ). Since the configuration of the stencil of any 
constant-coefficient scheme (like Eq. (2.1)) does not vary as j and n take different values, one 
concludes that the configurations of S 0 and So must be identical. In other words, B'A'D'C' can 
be made to coincide with ABCD by a simple translation in the x -direction (which, of course, is 
true). 

To further explore the significance of the above discussion, consider a constant-coefficient 
upwind scheme with the mesh point D missing from the stencil So (and thus D' missing from 
So). Under this circumstance, the configurations of So and So are different and it is obvious that 
the original equation and its image cannot satisfy the same upwind scheme. As a result, the 
scheme cannot be invariant under spatial reflection. 

Note that the spatial-reflection image of any mesh point (J,n), j,n =0, ±1, ±2, • • • , is (~j,n) 
(Here we assume that the plane of reflection is j = 0). Because the stencil So is formed by the 
spatial-reflection images of the mesh points of the stencil S 0 , So is the spatial-reflection image of 
S 0 . As a result of this observation, one can conclude from the above discussions that a constant- 
coefficient scheme cannot be invariant under spatial reflection if the configurations of its stencil 
and the spatial-reflection image are different. 

Similarly, one can conclude that a constant-coefficient scheme cannot be invariant under time 
reversal (space-time inversion) if the configurations of its stencil and the time-reversal (space- 
time-inversion) image are different. Note that the time-reversal (space-time-inversion) image of 
any mesh point (J,n ), j,n = 0, ±1, ±2, • • • , is (J, - n ) ((-;', - n )). 
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It has been shown that the Wendroff scheme is invariant under spatial reflection, time reversal, 
and space-time inversion. Another scheme which possesses the same properties is the Crank- 
Nicolson scheme [p.504, 3], i.e.. 


u" +1 - u] + -J [uj:l + u] +l - Uj-\ 1 ] = o 


( 2 . 11 ) 


Obviously the stencils of these two schemes satisfy the above necessary conditions for invariance 
under spatial reflection, time reversal, and space-time inversion. 


The Lax- Wendroff scheme [p. 101, 4], i.e., 

«+i _ v(v + l) 

Uj ~ 2 


«"-l + (1 -V 2 )mJ + 


v(v-l) 


M 7 +1 


( 2 . 12 ) 


is invariant under spatial reflection. However, its stencil does not meet the necessary condition 
for invariance under time reversal or space-time inversion. Generally, the stencil of an explicit 
upwind scheme does not meet any of the necessary conditions for invariance under spatial 
reflection, time reversal, and space-time inversion. 


3. VON NEUMANN ANALYSIS 

In this section, again we consider only the numerical analogues of Eq. (1.1). We assume that 
they are two-level linear difference schemes with real constant coefficients. As in Eqs. (2.1), 
(2.1 1), and (2.12), these coefficients are assumed to be functions of the Courant number v. 

For any one of the schemes referred to above, let G(v,0), ve /? ] , and < 0 < +°° be the 
amplification factor corresponding to the Fourier component e 1;0 . Here iav—l and 0 is the phase 
angle variation over Ax. Then, for any ve R \ and any 0, 

uj = [G(v,0)]V' e , j,n =0, ±1, ±2, • •• (3.1) 

is a solution to the scheme under consideration. Note that, given any v and 0, a two-level scheme 
has only one amplification factor. For any ve R i and any 0, let 

u] = [G(v,0)]V^ e , j,n = 0, ±1, ±2, • • • (3.2) 

be another solution to the same two-level scheme. Then one must have 

G(v,0) = G(v,0) , ve/? i, -~<0<+oo (3.3) 

By assumption, the coefficients of the scheme under consideration are real. It follows that the 
complex conjugate of a solution is also a solution. Thus Eq. (3.1) implies that, for any ve R i and 
any 0, 

uj = [G(v,0)]"e-‘ 70 , j,n =0, ±1, ±2, - (3.4) 

is also a solution to the same scheme. Since Eq. (3.4) will take the form of Eq. (3.2) if 0 is 
replaced by -0, Eq. (3.3) implies that 

G(v,-0) = G(v,0) , ve/?], — « < 0 < +°° (3.5) 

Let the scheme be invariant under spatial reflection, i.e., Eq. (2.8) is a solution to this scheme if 
Eq. (2.6) is. Then Eq. (3.1) implies that, for any ve R j and any 0, 

u] = [G(-v,0)]"e _iy0 , j,n = 0, ±1, ±2, • •• (3.6) 

is also a solution to this scheme. Since Eq. (3.6) will take the form of Eq. (3.2) if 0 is replaced by 
-0, Eq. (3.3) implies that 
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— oo < 8 < +oo 


(3.7) 


G(-v,-0) = G(v,0) , ve Ri , 
if the scheme is invariant under spatial reflection. 

Let the scheme be invariant under time reversal, i.e., Eq. (2.9) is a solution to this scheme if Eq. 
(2.6) is. Then Eq. (3.1) implies that, for any ve R i and any 0, 

uj = [(G(-v,Q)r'] n e ijd , j,n =0, ±1, ±2, (3.8) 

is also a solution to this scheme. By comparing Eq. (3.8) with Eq. (3.2) and using Eq. (3.3), one 
concludes that 

[G^v.e)]" 1 = G(v,0) , ve/?,, -oo< 8 <+oo (3.9) 

if the scheme is invariant under time reversal. 

Similarly, it can be shown that 

[GW.-e)]’ 1 = G(v,0) , ve/?! , -oo<0<+oo (3.10) 

if the scheme is invariant under space-time inversion. 

Note that any two of Eqs. (3.7), (3.9), and (3.10) implies the third. This is a result of the fact that 
any one of the three mappings SR, TR, and STI is the product of the other two mappings (in 
either order). 


By using Eqs. (3.5) and (3. 10), it is easy to show that 

G(v,0) = G(v,-0) = [ G(v,0) ] -1 (3.11) 

Thus we arrive at the conclusion that the scheme is neutrally stable, i.e., 

|G(v,0)| =1, ve/? i , -oo<8<+oo (3.12) 

if the scheme is invariant under space-time inversion. 


The amplification factors of the Wendroff scheme Eq. (2.1) and the Crank-Nicolson scheme Eq. 
(2.1 1), respectively, are 


G( _ cos(0/2) - i'vsin(Q/2) 
cos(0/2) + j'vsin(0/2) 


Wendroff 


(3.13) 


and 


1 - y sin0 

G(v,0) = Crank-Nicolson (3.14) 

iv 

1 + — sin0 
2 

For the Wendroff scheme, G(v,0) is not uniquely defined only when v = 0 and 0 = ±jc, ±3 re, 
±5rt, ■ ■ • . Thus we can choose R\ to be the set of all real numbers excluding 0. For the Crank- 
Nicolson scheme, G(v,0) is defined for all real v and real 0. Thus R i is the set of all real 
numbers. Since both schemes are invariant under spatial reflection, time reversal, and space-time 
inversion, both amplification factors satisfy Eqs. (3.5), (3.7), and (3.9) - (3.12). 


Contrarily, the amplification factor 

G(v,0) = 1 -v 2 (l - cos0 ) - /vsin0 -°°<v,0<+°° (3.15) 

of the Lax- Wendroff scheme Eq. (2.12) satisfies Eqs. (3.5) and (3.7) but not any one of Eqs. (3.9) 
— (3.12). This is consistent with the fact that the Lax- Wendroff scheme is invariant under spatial 
reflection but not under either time reversal or space-time inversion. 
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4. EXTENSIONS 


An extension to Eq. (1.1) is the convection-diffusion equation 

= 0 


du du d 2 u 
dt ° dx ~ * dx 2 


where a and p are arbitrary real constants. Eq. (4.1) is equivalent to 

du , du , d 2 u 

if 


dx 


dx 


> 2 


= 0 


(4.1) 


(4.2) 


(0 x' = -x , t' = t , a' = -a , and p' = p (4.3) 

or 

(«) x = x , a' = -a , and p' = -p (4.4) 

or 

(“0 x' = -x , t' = -t , a' = a , and p' = -p (4.5) 

In other words, Eq. (4.1) is mapped into itself under the mapping 

* -» x' , t -> t’ , a -> a' , and p — » p' (4.6) 

if x',t', a', and \i! are defined by any one of (i) — (Hi). Obviously, any one of the three mappings 
defined by Eqs. (4.6) and (4.3) - (4.5) can be considered as the product of the other two (in either 
order). With the same reasoning given in Section 1, Eq. (4.1) is said to possess invariant 
properties with respect to spatial reflection, time reversal, and space-time inversion. 


Let 


u = u 0 (x,t\a,ll) (4.7) 

be a solution to Eq. (4.1). By using the above invariant properties and a argument similar to that 
presented in Section 1 , it can be shown that 

u = u 0 (-x,t;-a,[L) (4.8) 

u = «o(x, -r;-a,-|x) (4.9) 

and 


are also solutions to Eq. (4.1). 


u = u 0 (-x,-t;a,-\i ) 


(4.10) 


Let 



2 pAr 
(Ax ) 2 


(4.11) 


Let u" be the dependent variables at the mesh point (J,n) of a constant-coefficient finite- 
difference analogue of Eq. (4.1). Let the mappings SR', TR\ and STI', respectively, be defined 
by (i) u] -> ulj (J,n = 0, ±1, ±2, • • • ) v -v, and y -> y, (u) u] -> u~ n (J,n = 0, ±1 , ±2, ■ • • ), 
v — > —v, and y — » — y, and (iii) u " — > (J,n =0, ±1 , ±2, • • • ), v — » v, and y — » — y. Obviously, 

any one of the mappings SR', TR', and STI' can be considered as the product of the other two (in 
either order). In this paper, an analogue of Eq. (4.1) is said to be invariant under spatial reflection 
(time reversal, space-time inversion) if the system of equations formed by this analogue is 
mapped into itself under the mapping SR' (TR', STI'). 
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Let R 2 be a set of ordered pairs of real numbers such that both (~x,y) and Qc, -y) are elements of 
R 2 if (x,y) is. For any ( v , y)g R 2 , let 

Uj = Ho(Jrny,y) , j,n=0, ±1, ±2, ••• (4.12) 

be a solution to a finite-difference analogue of Eq. (4.1). Then it can be shown that, for any 
(v,Y)€/? 2 > 


u] = «o(-J.«;-v,Y) . 

j,n= 0, ±1, ±2, • • • 

(4.13) 

u] = u Q (j, -n ;-v,-y) , 

j,n= 0, ±1, ±2, • 

(4.14) 

uj = moH,-«;v,-y) . 

j,n =0, ±1, ±2, • • 

(4.15) 


respectively, are solutions to this analogue if it is invariant under spatial reflection, time reversal, 
and space-time inversion, respectively. Note that the proof follows a line of arguments which 
was used to obtain similar results in Section 2. 

Also, it is obvious that the necessary conditions for invariance established in Section 2 regarding 
the stencil’s configuration remain valid for the current extension. 

Let the numerical analogue under consideration be a two-level linear difference scheme with real 
constant coefficients. These coefficients are assumed to be functions of v and y. Let G(v,y,0), 
(y,y)eR 2 and < 0 < +°°, be the amplification factor. Then by using a reasoning similar to 


that leading to Eq. (3.5), one concludes that 

G(v,Y,-0) = G(v,y, 0) , (v,Y)e Ri . < 0 < +°° (416) 

Similariy, by using Eqs. (4.13) - (4.15), it can be shown that (i) 

G(-v,y,-0) = G(v,y, 0) - (v,Y)e/?2. -°° < 0 < +°° (4.17) 

if the scheme is invariant under spatial reflection, (ii) 

[G(-v,-y,0)] _1 = G(v,y, 0) , (v,y) g /?2 » -°°<0<+°° (4.18) 

if the scheme is invariant under time reversal, and (iii) 

[G(v,-y,-0)] _1 = G(v,y, 0) , (v, Y)e/? 2 , -oo<e<+oo (4.19) 


if the scheme is invariant under space-time inversion. Obviously, any two of Eqs. (4.17) — (4.19) 
implies the third. This is a result of the fact that any one of the three mappings SR', TR', and 
STT is the product of the other two mappings (in either order). 

By using Eqs. (4.16) and (4.19), one has 

G(v,-y,0) = G(v,-y,- 0) = [G(v,y,0) ] -1 (4-20) 

As a result, one concludes that 

G(v,-Y,0) • G(v,Y,0) = 1 , (v,Y)eR 2 > -°°<0<+°° (4.21) 

if the scheme is invariant under space-time inversion. Note that Eq. (4.21) is reduced to Eq. 
(3.12) when y = 0. i.e., |X = 0. 

Two-level implicit schemes Eqs. (2.1) and (2.11) were given in Section 2 as examples of the 
numerical analogues of Eq. (1.1) which are invariant under spatial reflection, time reversal, and 
space-time inversion. Also we consider only two-level schemes in Section 3 and in the derivation 
of Eqs. (4.16) - (4.21). A three-level explicit scheme which is designed to solve Eq. (4.1) is the 
Leapfrog/DuFort-Frankel scheme [p.161, 4], i.e., 
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(1+Y)k > +1 = (l-Y )u] 1 + (V+Y )«"-1 - ( v — y) u]+i (4.22) 

This scheme is also invariant under spatial reflection, time reversal, and space-time inversion. It 
is reduced to (i) the Leapfrog scheme if y = 0 (i.e., |i = 0) and (ii) the DuFort-Frankel scheme if v 
= 0 (i.e., a = 0). 


The amplification factors of Eq. (4.22) are 

/- a. m 7cos0 - iv sin0 ± ^(ycosQ - ivsin8) 2 + 1 - -y 2 

G±(v,Y. 6) = L L_ t y*_i (423) 

Since (i) |cos8|<l and |sin0|<il, and (ii) G ± (O,O,0) = ±1, there exists a set R 2 of ordered pairs 
of real numbers such that (i) both (~x,y) and (x, -y) are elements of R 2 if (x.y) is, (ii) (x, -1) does 
not belong to R 2 for all x, and (iii) 

Re[G + (v,Y,8)] > 0 and Re[G_(v,Y,0)] < 0 (4.24) 

for all v, y, and 0 with (v,y)eR 2 and -« < 0 < + «,. Here Re[G ± (v,Y,0)] denotes the real part of 
G±(v,y f 0). Note that Eq. (4.24) is equivalent to 

Re[(G + (v,Y,0)r 1 ] >0 and Re[(G_(v,Y,0)r'] < 0 (4.25) 

Let (v,Y)e/? 2 - Let 

«" = [G(v, Y,0)]V' 9 (4.26) 

be a solution to Eq. (4.22). As a result of Eq. (4.24), one has 

f G+(v,y, 0) if Re[G(v,Y,0)] > 0 

G(v,Y. 9) = j _ (4.27) 

l G_(v,y, 0) if Re[G(v,Y,0)] < 0 

With the aid of Eqs. (4.24) — (4.27) and the fact that Eq. (4.22) is invariant under spatial 
reflection, time reversal, and space-time inversion, a line of arguments which were used to obtain 
Eqs. (3.5), (3.7), (3.9), (3.10), and (3.1 1) again can be invoked to show that 


G ± (v,y,- 8) = 

G±(v,Y.9) 

(4.28) 

G±(-v,y.- 8) = 

: G±(v,y, 8) 

(4.29) 

[G ± (-v,-y,0 )] -1 

= G ± (v,y, 8) 

(4.30) 

[G±(v,-y,-0 )] _1 

= G±(v,y, 8) 

(4.31) 

G±(v,-y, 8) • G ± 

(v,y. 0) = 1 

(4.32) 


for all v, y. and 0 with (v,y)eR 2 and -« < 0 < +oo. 

Note that 

G ± (v,O,0) = -ivsin0 ± ^1 - v 2 sin 2 0 (4.33) 

Thus | G ± (v,O,0) | =1 if | v | <1 . This result is consistent with Eq. (4.32). 

5. DISCUSSION AND CONCLUSION 

In this paper we study several invariant properties of the numerical analogues of Eq. (1.1). 
Particularly, it is shown that an arbitral two-level constant-coefficient finite-difference analogue 
of Eq. (1.1) is neutrally stable if it is invariant under space-time inversion. A similar study for 
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Eq. (4.1) is also presented. 


Since it is a common experience that the local behaviors of a nonlinear variable-mesh scheme 
may be predicted by using a local analysis (such as the von Neumann analysis) in which the 
dynamic coefficients and geometric parameters are frozen at their local values, die information 
gained and the techniques developed in the current study may also be useful in a similar study for 
a numerical analogue of nonlinear physical equations. 

As noted previously, the construction of a numerical scheme generally does not take into account 
the invariant properties of the physical equations to be solved. For an upwind scheme, the stencil 
contains more mesh points on the upwind side than on the downwind side. Thus the 
configuration of the stencil differs from those of its spatial-reflection image and space-time- 
inversion image. According to a discussion given in Section 2, this implies that the scheme is not 
invariant under either spatial reflection or space-time inversioa The stronger the upwind bias is, 
the further away this scheme tends to be from preserving these two invariances. According to 
analysis presented in Section 3, this also tends to increase numerical diffusion. 

For a two-level explicit scheme, the stencil generally contains several mesh points at the lower 
time level while only one at the upper time level. Thus the configuration of the stencil differs 
from those of its time-reversal image and space-time-inversion image. This implies that the 
scheme is not invariant under either time reversal or space-time inversion. Generally, the higher 
the order of accuracy of a scheme is, the more points will be at the lower time level of the stencil. 
In turn, this makes the scheme further away from preserving these two invariances. 

Currently, a new numerical framework for solving conservation laws — the method of space- 
time conservation element and solution element is being developed by Chang and To [2] . This 
framework is fundamentally different from the well established methods, i.e., finite difference, 
finite volume, finite element, and spectral methods, in both concept and methodology. It may be 
used to solve both inviscid and viscous flow problems. 


A two-level explicit numerical analogue of Eq. (4.1) was constructed using this framework [2]. 
An unique feature of this scheme is that there are two dependent variables and two equations 
associated with each solution element. These two variables may be considered as the numerical 
analogues of u and du/dx. Because of this feature, contrary to traditional two-level explicit 
schemes, it is invariant under spatial reflection, time reversal, and space-time inversion (to be 
shown in another paper). Its two amplification factors are [o± (v,y, 0)] with 


o ± (v,y,9) ^ 


ycos^) - i'vsin(-|-) ± [ y cos( — ) - tvsin(^)] 


+ 1 - Y 


(5.1) 


1 + Y 


where Y 2y/(1-v 2 ) and 1-v 2 > 0 [2]. A comparison between Eqs. (5.1) and (4.23) reveals 
that the expression on the right side of Eq. (4.23) can be converted to that on the right side of Eq. 
(5.1) if Y. v, and 0, respectively, are replaced by y> v, an( * 0/2- Because of this remarkable 
similarity, it can be shown [2] that the stability condition of the new scheme, as in the case of the 
Leapfrog/DuFort-Frankel scheme, is essentially the CFL condition and thus independent of the 
viscosity coefficient |X. Therefore, the new scheme is unconditionally stable in the case of pure 
diffusion. Also, as in the case of the Leapfrog/DuFort-Frankel scheme, the new scheme has no 
numerical diffusion in the absence of viscosity. 
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